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Abstract. We present the results of two simultaneous X-ray/radio observations of the low-mass X-ray binary GX 
13+1, performed in July/ August 1999 with the Rossi X-ray Timing Explorer and the Very Large Array. In X-rays 
the source was observed in two distinct spectral states; a soft state, which had a corresponding 6 cm flux density 
of ~0.25 mjy, and a hard state, which was much brighter at 1.3-7.2 mjy. For the radio bright observation we 
measured a delay between changes in the X-ray spectral hardness and the radio brightness of ~40 minutes, similar 
to what has been found in the micro-quasar GRS 1915+105. We compare our results with those of GRS 1915+105 
and the atoll/Z-type neutron star X-ray binaries. Although it has some properties that do not match with either 
atoll or Z sources, GX 13+1 seems more similar to the Z sources. 
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1. Introduction 

Based on their correlated spectral and variability prop- 
erties, the brightest persistent neutron star low-mass X- 
ray binaries (LMXBs) ar e often divided in two groups: 
the atoll and Z so urces ijHasinger fe van der Klisl il989: 
van der Klislfl995^ . after the tracks they trace out in 
X-ray colour-colour diagrams (CDs). Although atoll and 
Z sources share some variability and spectral proper- 
ties (in the X-ray band), there are significant differ- 
ences between the two groups: atoll sources are less lu- 
minous, have harder X-ray spectra and s how stronger 
rapid time variability than the Z sources l|van der Klisl 
Il995al). Atoll sources ar e also less luminous in the radio 
llFender fc Hendrvll2000lh Although some of the observa- 
tional differences can be accounted for by differences in 
the mass accretion rate, with the Z sources probably ac- 
creting near the Eddington rate and atoll sources at rates 
ranging from near Eddington to less than 10% of it, it is 
generally believed that additional differences, e.g. in the 
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neutron star proper ties, are required to explain all obser- 
vational differences l|Hasinger fc van der KlislflOS^ . 

The nature of GX13+1, one of the brightest neu- 
tron star LMXBs, is still ambiguous. Schulz et al. I <ll989f) 
group ed it with the high luminosity sources i|Schul^et^d1 
1989), which included the six sources that were later la- 
beled as Z sources. IHasinger fc van der Klisl lll989h classi- 
fied it as a bright atoll source, although they noted that 
of all atoll sources it showed properties which were clos- 
est to those seen in the Z sources, most notably its flar- 
ing branch-like appearance in the CD and the featureless 
power law noise in the power spectrum that is typical for 
that spectra l state. Although the EXOSAT observations 
analyzed bv lHasinger fc van der Klisl l)1989j) only showed 
the source in the so-called 'banana branch ' state, the 
bimod al behavior of the source reported by IStella et alJ 
(1985) strongly suggests that the source occasionally en- 
ters a different state. The latter seemed to be confirmed 
by the first observations of the source with the Rossi 
X-ray Timing Explorer (RXTE), in which a clear two 
branc hed structure was found in the CD l|Homan et alJ 
1998). Although the pattern in the CD resembled both Z 
and atoll source tracks, the resemblance to Z sources was 
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strengthened by the discovery (in the same observations) 
of a 57-69 Hz quasi-periodic oscillation (QPO), which had 
properties similar to that of the horizontal branch QPO 
in the Z sources. The CD of GX 13+1 in a recent paper 
bv lMuno et alJ l|200lj) . which for the first time displays a 
'complete' pattern, shows a sharp vertex, similar to the 
normal branch/flaring branch vertex seen in the Z sources 
and quite unlike the ra ther smooth curves seen in atoll 
sources. More recently, ISchnerr et alJ l|2003h analyzed a 
large set of RXTE data and conclude that many of the 
source's properties do not fit within the atoll/Z frame- 
work, although they favor the option of an atoll source. 
They suggest that part of its unusual behavior can be 
explained with the presence of a relativistic jet, that is 
almost pointing directly towards us. 

Observations in the infrared suggest the pres- 
ence of a K giant seco ndary I Garcia et al.1 Il992l 
iBandvonadhvav et al . 1999) and a possible orbital or 
precessional modulation wi th a period of —20 days 
l|Bandvopadhvav et al.l Eo02j) ; these are properties that 
are thought to be more typical for Z sources than for 
atoll sources. G X13+1 also shares a common mean ra- 
dio l uminosity l)Grindlav fc Seaciuistl Il98f3: lOarcia et alJ 
Il988l) with the Z sour ces and black hole candidates 
Render fc Hendrvl |2000() . suggesting a similar origin for 
the quiescent radio emission from persistent black hole and 
Z source X-ray binaries, and GX 13+1. None of the (other) 
atoll sources is consistent with this relation. Although the 
source is variable in radio on time scales of less than an 
hour, no clear r elation between X- ray and radio luminosi- 
ties was found l|Garcia et al.lll988j) . 

In this paper we present the results of a coordinated 
X-ray/radio campaign, which had as its principal aim to 
investigate the possibility that GX 13+1 shows a similar 
X-ray state dependence of its radio emission as is found 
in the Z sources. In most of those sources the radio lumi- 
nosity decreases from the horizontal branch to the flaring 
branch. We find similar behavior in GX 13+1. 



2. Observations and analysis 

2.1. X-ray observations 

Our X-ray data were obtained simultaneou sly with the 
the P roportional Counte r Array (PCA; IZhang et alJ 
Il993l I.Tahoda et all Il996l) and th e High Energy X- 



ray Timing Experimen t (HEXTE; iGruber et all 11996! 
Rothschild et al. | [l998) onboard RXTE llBradt et al 



1993). The observations were performed in 1999 between 



July 31 23:52 UTC and August 1 10:39 UTC (obs. 1), and 
on August 04 between 02:44 and 12:05 UTC (obs. 2). The 
total exposure times for the two observations were, respec- 
tively, -22.5 ks and -18.3 ks. The PCA and HEXTE data 
were obtained in several different modes; the spectral and 
timing properties of these modes are given in Table ^ For 
all modes we discarded data taken during Earth occulta- 
tions and passages through the South Atlantic Anomaly. 



The Standard 2 data were used to produce light 
curves, colour curves, a colour-colour diagram (CD), and a 
hardness-intensity diagram (HID), and to perform a spec- 
tral analysis. Only data from Proportional Counter Units 
(PCUs) and 2 were used (all layers), since these were 
the only two that were active during all our observations. 
The data were background subtracted; dead time correc- 
tions (—3-3.5%) were only applied for the spectral anal- 
ysis. A soft colour was defined as the ratio of count rates 
in the 4.2-7.5 keV and 2.5-4.2 keV energy bands, and a 
hard colour as the ratio of count rates in the 10.0-18.5 
keV and 7.5-10.0 keV energy bands (these four energy 
bands correspond, respectively, to Standard 2 channels 
(running from 1 to 129) 7-14, 3-6, 21-40, and 15-20). A 
CD and a HID were produced by, respectively, plotting 
hard colour versus soft colour and hard colour versus the 
2.5-18.5 keV count rate, for each 16 s data point. The 
soft and hard colour curves, which had an initial time 
resolution of 16 s, were rebinned to a time resolution of 
1024 s, to allow a better study of their variations on long 
time scales. The PCA spectra were created using the stan- 
dard FTOOLS V5.2 routines. Systematic errors of 0.6% 
were added and response matrices were produced using 
PCARSP (V8.0). The details of the models used to fit the 
spectra of GX13+1 are discussed in Section ITTI 

HEXTE light curves were produced by running 
the FTOOLS V5.2 script hxtlcurv (which automat- 
ically corrects for background and deadtime) on the 
E_8us_256_DXlF mode data of cluster A. Additional rebin- 
ning was applied, resulting in 1024 s time bins, to achieve 
better statistics. 

The rapid X-ray time variability of the source was 
studied in terms of power spectra. For this analysis we 
performed fast fourier transforms (FFTs) of the high time 
resolution (1/8192 s) data . Power spectra were produced 
with frequencies of 1/1024-512 Hz (2-26.5 keV band) and 
1/16-2048 Hz (2-26.5 keV, 2-6.3 keV, and 6.3-26.5 keV 
bands). No background or dead time corrections were ap- 
plied to the data prior to the FFTs; the effects of dead time 
were accounted for by our fitting method. The power spec- 
tra were selected on time, count rate and/or colour and 
subseque ntly rms normalized a ccording to a procedure de- 
scribed in Ivan der Klisllll995bh . The resulting power spec- 
tra were fitted with a constant, to account for the dead- 
time modified Poisson level, a power law, a zero-centered 



Table 1. Data modes for the RXTE/PCA and 
RXTE/HEXTE. 



Mode 


Time 


Energy 


Energy 




res. (s) 


range (keV) 


channels 


Standard 1 


1/8 


2-60 


1 


Standard 2 


16 


2-60 


129 


SB_125us_0_13_ls 


1/8192 


2-5.9 


1 


SB_125us_14_17_ls 


1/8192 


5.9-7.6 


1 


E_125us_64M_18_ls 


1/8192 


7.6-60 


64 


E_8us_256_DXlF 


1/131072 


10-250 


256 
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Fig. 1. The evolution of GX 13+1 during the observations on 1999 August 1 (left column) and 1999 August 4 (right 
column). a,b: PCA Count rate in the 2-25 keV range (16 s bins). c,d: HEXTE 20-100 keV count rate (rebinned to 
1024 s). e,f: Soft color (rebinned to 1024 s). g,h: Hard color (rebinned to 1024 s). i,j: VLA 6 cm radio flux density 
(1800-4800 s [i] and 900 s [j]). Errors on the flux density in panel i are of the order of 50 /iJy. For the PCA count 
rate and colours only data from PCUs and 2 were used. For the HEXTE count rate only Cluster A was used. Note 
that all quantities have the same ranges for obs. 1 and 2, except for the hard colour (to allow an easier comparison 
the with changes in the radio flux density). 



Lorcntzian, to account for the deviations from the power 
law noise, and a Lorentzian for a weak QPO in the second 
observation. Errors on the fit parameters were determined 
using Ax 2 = 1. Upper limits on QPOs were determined 
by fixing the frequency and FHWM to (a range of) values 
and using A\ 2 = 2.71 (95% confidence). 



2.2. Radio observations 

The radio data were obtained with the Very Large Array 
(VLA) radio observatory in its most extended (A) config- 
uration.. GX13+1 was observed at 6 cm on 1999 August 
1 between 00:33 and 09:31 UTC (obs. 1), and on 1999 
August 4 between 00:22 and 08:51 UTC (obs. 2). During 
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Fig. 2. (a) Colour-colour diagram and (b) hardness- 
intensity diagram for the two RXTE/PCA observations 
of GX 13+1. Points of the first observations are depicted 
by the open circles, those of the second observation by 
crosses. Each data point represent a 16 s interval. Typical 
error bars are shown in the lower-right corners. 



the first observation 25 of the 27 VLA antennas were used; 
26 were used during the second observation. Due to tech- 
nical problems ^39% (obs. 1) and ~8% (obs. 2) of the 
total observing time was lost. Observations were made 
simultaneously in both circular polarizations in each of 
two independent 50 MHz bands centered on 4885.1 and 
4835.1 MHz. Flux densities, which are all Stokes I, were 
referenced to those of 1328+307 (3C 286), taken to be 
7.462 and 7.510 Jy at 4885.1 and 4835.1 MHz respectively 
(Perley et al., priv. comm.). The overall flux density scale 
is probably good to at least ~ 5%. To calibrate the com- 
plex antenna gains the standard calibration source 1817- 
254 was observed at 30 minute (obs. 1) or 6 minute (obs. 
2) intervals; the reason for the denser sampling of the cali- 
bration source in obs. 2 was that during obs. 1 a large part 
of the data was rendered useless due to a combination of 
bad weather conditions and too long intervals between 
calibration observations. The data were reduced and an- 
alyzed using the Astronomical Image Processing System 
(AIPS). For each day, images were first made using the en- 
tire data sets, ignoring any variability of GX 13+1. These 
images were used to find nearby confusing sources, which 
were then subtracted from the original uv-data. These uv- 
data were then split into 15-minute bins (or longer, in the 
case of weak signals) and imaged to give the flux density 
history of GX 13+1. This procedure allows the full synthe- 
sis mapping of confusing sources, while retaining the best 
possible sensitivity to fluctuations of GX 13+1. The flux 
densities are the average of the of maximum flux density 
and the integrated flux density in the region containing 
the source in the cleaned image. Note that because of cal- 
ibration and other "dead" time, the amount of data in 
the each bin is usually less than the sampling time. Errors 
on the flux density were conservatively calculated as the 
sum of the squares of (1) the difference of the maximum 
and the integrated flux density and (2) the measured off- 
source rms. Early in obs. 1 and during the beginning and 
end of obs. 2 the source wandered away from the image 
center, probably because the elevations of GX 13+1 and 
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Fig. 3. The ratio of the mean raw RXTE count rate spec- 
trum of observation 2 to that of observation 1, revealing an 
overall increase and a relative hardening of the spectrum. 

the calibrator were rather different there - the effects of 
these problems can clearly be seen in the increase of the 
error bars at the end of the radio light cure in Fig. [T], 

3. Results 

3.1. X-ray observations 

From Figs. it is clear that the X-ray properties of the 
source changed between the first to the second observa- 
tion. The second observation has a higher count rate, both 
in the 2-25 keV (Fig. [TJi and b) and 20-100 keV bands 
(Fig. and d) and shows more variability in the 2-25 
keV band. Also, the second observation tends to be spec- 
trally harder; while the PCA count rates increased only by 
a factor of ~1.3, the HEXTE count rates increased by a 
factor of ^3.2. The spectral difference is most clearly seen 
in the CD and HID shown in Fig. |3 where the two ob- 
servations show up as distinct patches, and Fig. |3 which 
shows the ratio of the spectra of obs. 2 and 1. Note that 
the small bridge between the two patches in the HID (Fig. 
2b) does not represent a real connection between the two 
patches, but corresponds to the upper right part of patch 
traced out during obs. 1 in the CD. It is not clear from 
our observations how the source bri dges the gap between 
the two patches in HID - we refer to lSchnerr et al.l ((2003) 
for an observation that shows a transition between the 
two patches. Although type I X-ray bursts have been ob- 
serve d from GX 13+1 l|Fleischmanlll985l iMatsuba et alJ 
1995), none were seen during our observations. 

3.1.1. X-ray spectra 

X-ray spectral fits were performed with XSPEC l|Arnaud| 
Il99fl V11.2) to the 3-25 keV PCA spectra using two dif- 
ferent models. The first model i s base d on the contin- 
uum model used by llJeda et all ll200ll). for their ASCA 
1-10 keV spectra of GX13+1, and ISidoli et~aH ll2002ft . 
for their 2-10 keV XMM-Newton spectra. It consists of 
a black body (bbody in XSPEC), and a disk black body 
(diskbb). Galactic neutral hydrogen absorpt ion (phabs) 
was fixed at Nh = 2.9 10 22 atoms cm~ 2 (|Ueda et alJ 
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Table 2. Spectral fit results. Errors on the fit parame- 
ters are 68% confidence limits, upper limits represent 95% 
confidence limits. The 3-10 and 10-25 keV fluxes are un- 
absorbed. 



Parameter 


Obs. 1 


Obs. 2 


CN 


Njj (atoms cm ^) 


z.y iu ^nxea ) 


z.y iu ^nxea j 


C 
D 
CD 


KLw ^e v ) 


n oQ-i-A no 
U.yozhu.Uz 


c\ oo-i-fi m 
U.yzzhU.Ul 


E 
\ 






04-1-0 04 


V) 

E n 


T 


7.2±0.2 


8.10±0.15 


^ 


Rw (km) 


21.7±0.9 


27.2±0.7 


* 
>> 


Efe (keV) 


6.42 (fixed) 


6.42 (fixed) 




rz 
a) 


FWHM (keV) 


0.78±0.13 


4 +0 " 2 


O" 

(D 


Fe EW (eV) 


223 


<96 


Li_ 


E edge (keV) 


9.1±0.1 


9.1 (fixed) 






0.12±0.03 


<0.03 


O 


-F3-10 (ergs cm^s" 1 ) 


1.5 10" 8 


2.0 




fio-25 (ergs cm _2 s _1 ) 


1.8 1(T 9 


3.6 10" 9 




xled (dof) 


0.6 (40) 


0.9 (41) 





Fit parameters of the comptt model - kTw: input soft photon 
(Wien) temperature - kT e : plasma temperature - t: plasma 
optical d epth - Rw'- effective Wien radius for the soft seed 
photons qin 't Zand et al.lll999l ). here derived for an assumed 
distance of 7 kpc 



120011) . A reasonable fit was found for obs. 2 (Xred = 1-52 - 
with an unphysically small disk radius), but not for obs. 1 
(.Xred — H-3). Adding a line (gauss) arou nd 6.4 keV 
and a n absorption edge (edge) around 8 keV l|Ueda et alJ 
l200lh did not lead to acceptable fits. The second model 
we tried was a t hermal Comptonization model (comptt, 
iTitarchukl Il994l iHua fe Titarehukl Il995|) . which has re- 
cently been applied successfully to the broad band con- 
tinuum of several bright neutron sta r LMXBs (see, e.g., 
Idi Salvo et a"fll2000l llaria et al.ll200lh . The N H was again 
fixed to Nh = 2.9 10 22 atoms cm' 2 . For the second ob- 
servation a reasonable fit was obtained (xl e d = ^-^)^ but 
the first observation showed large residuals between 5 and 
10 keV (x 2 ed =15.6). Adding a line around 6.4 keV and an 
edge around 9 keV greatly improved the fit (x 2 ed =0.6). 
For the second observation the addition of the line and 
edge led to a small improvement (Xred = ^-^)' with the two 
components not being detected significantly. The best-fit 
results with the four-component model are given in Table 
El Some caution should be taken with interpreting the 
spectral results, as the low line energy (which tended to 
decrease when not fixed) and the high edge energy indicate 
different ionization stages of Fe. Also the large spectral 
changes within obs. 1 may have been partly responsible 
for the observed residuals. For both observations no black 
body component was needed, as the Comptonization com- 
ponent provided good fits at low energies. 

For a more model independent comparison of the spec- 
tra of the two observations we plot the ratio of observa- 
tions 2 and 1 in Fig. [3J While below 7 keV the ratio is 
rather constant, above that energy the spectrum of obs. 2 
becomes increasingly hard. 




0.1 1 10 

Frequency (Hz) 

Fig. 4. Power spectra (2-26.5 keV) of the two 
RXTE/PCA observations (lower: obs. 1, upper: obs. 2) 
in an vP(v) representation. Note that the upper power 
spectrum has been multiplied by a factor of 10, to avoid 
overlapping of the power spectra. The Poisson level has in 
both cases been subtracted for plotting purposes. 



3.1.2. X-ray variability 

From the PCA light curves in Fig.QJit is clear that, at least 
on time scales of minutes to hours, the source was more 
variable during the second observation. The 2-26.5 keV 
power spectra of the two observations are shown in Fig. 
01 the results of the power spectral fits are given in Table 
131 Both power spectra are relatively featureless and their 
continuum can be well fit by a combination of power law 
noise and a zero-centered Lorentzian. While the power law 
noise has an identical index in the two observations, the 
FHWM of the two zero-centered Lorentzians, and hence 
their maxima in a vP v plot, differ by more than two orders 
of magnitude, suggesting that they might not be related. 
We added a Lorentzian to the fit function of the second 
observation, for a QPO at 60 Hz - this QPO has a sta- 
tistical significance of 3.2er. No similar QPO was found in 
the first observation. We also measured the total rms in 
two frequency intervals; 0.001-1 Hz and 1-100 Hz. In both 
frequency intervals the fractional variability is stronger in 
the second observations, with the difference being more 
pronounced at high frequencies. We searched for high fre- 
quency QPOs in the 1/16-2048 Hz power spectra with 
the FWHM fixed to 150 Hz, but only upper limits could 
be determined. For obs. 1 they were: 2.1% (2-26.5 keV), 
2.8%(2-6.3 keV), and 4.0% (6.3-26.5 keV) and those for 
obs. 2: 2.2% (2-26.5 keV) , 3.7% (2-6.3 keV), and 3.6% 
(6.3-26.5 keV). 
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Table 3. Power spectral fit results. As a model indepen- 
dent measure of the variability we also include the total 
fractional rms in the 0.001-1 Hz and 1-100 Hz ranges. 
FWHM stands for full-width-at-half-maximum. 



Parameter 


Obs. 1 


Obs. 2 


Power law rms a (%) 


5.15±0.12 


7.71±0.11 


Power law index 


1.09±0.01 


1.09±0.01 


Lorentzian rms 6 (%) 


2.99±0.01 


3.4±0.6 


Lorentzian FWHM (%) 


0.45±0.04 


100±30 


Lorentzian freq. (Hz) 


(fixed) 


(fixed) 


QPO rms c (%) 


<1 


2.1±0.4 (3.2o-) d 


QPO FHWM (Hz) 


18 (fixed) 


18±6 


QPO frequency (Hz) 


60 (fixed) 


60±2 


rms [0.001-1 Hz] (%) 


6.6±0.02 


9.15±0.02 


rms [1-100 Hz] (%) 


3.17±0.10 


5.93±0.06 


Xled(d.0.f) 


1.07 (209) 


1.01 (210) 



a Integrated between 0.001 and 1 Hz b Integrated between 
and +oo Hz c Integrated between — oo and +oo Hz 
Significance is calculated from the power, not the fractional 
r.m.s. 



3.2. Radio observations 

During the first observation the source was only sometimes 
detected significantly, at a very low level and with low 
variability; the average flux density was ~0.25 mJy. The 
radio flux density in the second observations was much 
higher, varying between 1.3 and 7.2 mJy. The flux density 
showed a slow increase for about 4.5 hours, followed by 
a gradual decrease. No significant variability was detected 
on time scales less than 15 minutes. An upper limit on the 
size of the source of 0.15 arcsecond was obtained. 

3.3. X-ray/radio comparison 

Comparing panels i and j of Fig. ^with panels a and b of 
the same figure shows that, apart from the fact that the 
second observation is brighter both in X-rays and radio, 
there is no apparent relation between X-ray count rate 
and radio flux, at least not on the time scales on which the 
latter varies. On the other hand, a comparison of the radio 
flux and hard X-ray colour of obs. 2 reveals similar broad 
peaks centered around ^05:00 UTC. As is obvious from 
panels nji and^ there exists a significant lag between the 
peaks of the hard colour curve and the radio light curve. 
Fitting the broad peaks in the hard color (02:52-06:02 
UTC) and radio curves (03:53-06:53 UTC) both with a 
Lorentzian gives a lag between the two peaks of 42±3 
minutes. 

4. Discussion 

We have observed the neutron star LMXB GX 13+1 simul- 
taneously in X-rays and radio. The source was found in 
two clearly distinct X-ray states; the spectrally hard state 
had associated radio fluxes that were between 4 and 18 
time higher than the maximum we detected in the softer 



X-ray state. This dependence of the radio flux on X-ray 
state is similar to what is found for other bright neutron 
star LMXBs (the Z sources; see iHiellming Rr. Ha.nl H 995l 
and references therein) and more re cently also in the less - 
luminous atoll source 4U 1728-34 l|Migliari et al.ll2003h . 
More specifically, our observations strongly suggest that 
the radio flux of the source is related to the X-ray spectral 
hardness, on time scales of both hours and days. 

4.1. Outflow 

It is generally believed that the radio emission from neu- 
tron star and black hole X-ray binaries is produced by 
highly relativistic electrons that interact with magnetic 
fields to produce synchrotron radiation. High resolution 
radio observations of a handful of X-ray binaries show 
that these electrons reside in powerful, collimated out- 
flows, commonly referred to as jets. It is assumed that 
in the X-ray binaries for which jets have not (yet) been 
directly observed the radio emission originates in a similar 
outflow. To see whether this could also be the case for GX 
13+1, we estimate the size of the emission region. If, for 
an assumed spherical region, this size is larger than the bi- 
nary separation, the highly relativistic synchrotron plasma 
cannot be contained in the system and the most plausible 
option would then be that we are dealing with an out- 
flow. Assuming a maximum brightness temperature T r of 
< 10 12 K (see e.g. lKellermann fc Paulinv-Totldll969T) and 
a distanc e of 7 kpc llBandvopadhvav et al.ll999(l we derive, 
following iFenderl l|2003f) . a minimum size for the emitting 
region of ^15 Rq (10 12 cm) at the peak of the radio flare 
in obs. 2. A size estimate based on the fastest part (1.5 
hr) of the rise of the radio flare (< cAt ~ 1.6 ■ 10 14 cm) is 
consistent with this value. Even for a system with equal 
masses for the primary and secondary and an orbital pe- 
riod of 20 days the emission region is at least comparable 
and probably larger than the size of the binary system 
(~ 3 • 10 12 cm), suggesting an outflow is also present in 
GX 13+1. 

4.2. X-ray/radio connection 

Although the exact mechanism for producing jet outflows 
in X-ray binaries is still not clear, the energy needed to 
achieve the inferred observed high bulk velocities suggests 
that they form in the inner parts of the accretion flow, 
where also most of the X-rays are produced. Simultaneous 
X-ray and radio observations of GRS 1915+105 seem to 
confirm such a dir ect link between the outflow and in- 
ner a ccretion disk ((Mirabel et al.lll998l: iKlein-Wolt et alJ 
20021). Clear patterns of X-ray/radio behavior are also ob- 
served in other black hole and Z source X-ray binaries; 
in general the spectrally hard X-ray states have a higher 
radio luminosity than the spectrally soft states. 

Our observations of GX 13+1 reveal a similar pattern, 
with the second observation being both more luminous in 
the radio and showing a harder energy spectrum (Figs.^ - 
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Table 4. Observed mean and maximum radio flux den- 
sities and estimated distances for the six Z sources and 
G X 13+1. This table is rep roduced from Tables 2 and 3 
in l|Fender fc Hendrvll2000ji . with updated values for GX 
13+1. 



Source 


Mean a 
(mJy) 


Max 6 
(mJy) 


distance 
(kpc) 


Refs 


Sco X-l 


10 ±3 


22 


2.0+1.0 


1,2 


GX 17+2 


1.0 + 0.3 


13.4 


7.5 + 2.3 


1.3 


GX 349+2 


0.6 + 0.3 


1.3 


5.0 + 1.5 


4,5 


Cyg X-2 


0.6 + 0.2 


3.4 


8.0 + 2.4 


1,6,7 


GX 5-1 


1.3 + 0.3 


1.6 


9.2 + 2.7 


1 


GX 340+0 


0.6 + 0.3 


0.6 


11.0 + 3.3 


8 


GX 13+1 


1.8 + 0.3 


7.7 


7+1 


9-13 



a Mean cm radio flux density 

b Maximum radio flux density at 6 cm 



Refs 1: IPenninxl (ll98Sf) . 2: iBradshaw et all Jl99St) 3: 
Penninx et alJ dl98gl). 4 : ICooke fc PonmarJ 199ll ). 5: 
ClirisTlan^^^wankl ll997l) . 6: IHiellming et ali 



7: ICowlev et al] J197L.. 
Grindlav fc Seaauistl fl986) 



8: |Pejminx^^L| 

10: iGarcia et al] 
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Berendsen et al.1 J2000T) . 12: iBandvopadhvav et al.1 J1999T) . 13 

this work 



While the change between the soft and hard spectral 
states of GX 13+1 most likely occurred on a time scale 
of one or two days, we also find a relation between the 
X-ray and radio properties on a time scale of a few hours: 
a hardening of the X-ray spectrum was followed by an 
increase in the radio luminosity with a delay of ~40 min- 
utes. This is the first time such a short-term X-ray/radio 
connection and delay have been found and measured in 
a neutron star LMXB. It is interesting to note that a 
similar delay has also been found on several occasions 
in the galactic black ho le X-ray binary GRS 1915+105 
l|Klein-Wolt et alJ [2002) . In that case however, the de- 
lay (~45-60 minutes) was with respect to the start of 
the dips in the RXTE/PCA count rate. These dips co- 
incided with a strong spectral hardening. Also, the X- 
ray dips and radio flares in GRS 1915+105 were a rapid 
recurring phenomenon, whereas the r adio flare in GX 
13+1 was more likely a singular event. iKlein-Wolt et all 
l)2002f) explained the observed delay in GRS 1915+105 as 
the time it takes for the flow to become optically thin 
in the radio; this eff ect is clearly observed in Fig. 9 of 
Dha wan et al . (2000), which shows light curves of the di- 
rectly imaged jet of GRS 1915+105 at different radio 
wavelen gths that peak later w ith increasing wavelength 
sec alsolPoolev fc Fenderlll997|) . Observations of Sco X-l 



( Fomalo nt et alJl2001l) suggest that such a delay can also 
result from the transfer time of the energy from the core 
to the radio lobes. 

The ratio of (radio flux density change) /(X-ray flux 
change) is much larger in GX 13+1 than in GRS 
1915+105. If we assume that the inner disk is ejected to 
form the outflow, as has been proposed for GRS 1915+105 



to explain the X-ray variations, our second observation of 
GX 13+1 suggests that a smaller fraction of the inner disk 
is ejected in this case. The fact that clear dips are ob- 
served in GRS 1915+105 might be related to the different 
nature of the compact object, or to a much l arger amount 
of matter being expelled. Following iFenderl {2003) we can 
estimate the mean power of the radio event in obs. 2 to 
be ~ 1.8 • 10 36 erg s _1 . While this is only ~1 percent of 
the 3-25 keV X-ray luminosity (~ 1.4 ■ 10 38 erg s" 1 ), the 
kinetic energy associated to the (possibly relativistic) out- 
flow may increase this number to a larger fraction of the 
total accretion energy. 

Although short term variations in the radio luminosity 
of GX 13+1 have been found bef ore, they were only com- 
pared to the X-ray count rates ijGarcia et al.lll98q) . and 
not to the X-ray spectral properties. No correlations be- 
tween radio luminosity and X-ray count rate were found, 
most likely because X-ray data were only available for the 
radio weak part of their data se t. The time scale of s hortest 
radio fluctuations observed bv lGarcia et"afl l|l988|) (a few 
hours) is consistent with that of the flare in our second ob- 
servation. The 6 cm peak flux density they measured was 
2.2 mJy, about a factor of ~3.3 lower than our maximum 
flux density. Interestingly, when fitting our X-r ay spectra 
with the same Nh used bv lGarcia et alJ <jl988[) and using 
their assumed distance of 7 kpc we find a similar difference 
for the 1-20 keV luminosity. 



4.3. Comparison with Z and atoll source observations 

Before discussing the nature of GX 13+1 in view of the 
atoll/Z source classification we compare its X-ray/radio 
behavior with that o f sources from both classes. Although 
ISchnerr et aD l)2003|) concluded that GX 13+1 is neither 
a Z nor an atoll source, one could, purely based on its ap- 
pearance in the CD, argue that the source was in the atoll 
source banana branch and island state, during obs. 1 and 
2 respectively, or the Z source flaring branch and normal 
branch. 

Bas ed on observing cam paigns of th e Z sources GX 
17+2 l|Penninx et all ll98Sl). Cvg X-2 faiellming et alJ 
Il990al). Sco X-l llHiellming et alJ ll990rJ) and GX 5-1 
l|Tan et alJll992j) . IPenninxl <|l989|7 suggested that all Z 
sources share a common radio, UV, and X-ray lumi- 
nosity on thei r normal branch. Subse quent detections 
of GX 349+2 f Cooke & PonmarJ Il99l and GX 340+0 
l|Penninx et al.lfl993|) at approximately t he predicted ra- 
dio brightness seemed to confirm this idea. lPennimd |l9891 
derived a normal branch luminosity of ~ 1.6 • 10 38 erg s _1 
in the 1.5-15 keV band. In the same energy band we ob- 
tain (extrapolating our fit to lower energies) ~ 1.7 • 10 38 
erg s" 1 , which is remarkably close the value of IPenninxl 
lll989h . 

The assumption of lPenninxl l|l989jl might be not com- 
pletely valid; u sing the distance es t imate s and average flux 
densities from iFender fc Hendrvl {2000} (see Table ^ it 
seems that there is a considerable spread in the average 
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radio luminosity of the six Z sources. Moreover, none of 
them has a higher value than GX 13+1. However, this 
average radio luminosity depends strongly on the time a 
source spends in a radio bright state. Taking the maxi- 
mum radio flux densities (at 6 cm) reported in the liter- 
ature (see references in Table 0J should partly correct for 
this, which results in GX 17+2 being the most luminous 
with GX 13+1 being second. The overall behavior of the 
radio brightness along the track in the CD is also similar 
to that observed in several of th e Z sources l|Penninx et alJ 
ll988llHieilming et al.lll990albl) . 

The only atoll source for which a clear pattern in 
the X-ray /radio em ission has been found is 4U 1728-34 
( Mi gliari et alJl2003j) . In that source the radio flux density 
was highest in the island state; it increased by a factor ~6 
from the hard part of the island state toward the transi- 
tion between the island state and, what seemed to be, the 
(softer) banana branch. After this transition it dropped 
sharply by a factor of ~6. While the decrease of radio flux 
with spectral hardness (in the island state) is opposite to 
overall behavior in GX 13+1 and several Z sources (where 
we see an increase with spectral hardness along the track 
in the CD), the factor of ~6 difference between the island 
state and banana branch is similar in sign and magnitude 
to the difference between our second and first observation. 
The average X-ray luminosity of 4U 1728 -34 (assuming a 
distance of 5.2 kpc I Galloway et alJEo02l) ) is more than a 
factor 10 lower than that of GX 13+1, whereas the 6 cm 
flux density is probably more than a factor 100 lower. 

4.4. Z or atoll? 

The nature of GX 13+1 and its place in the Z/atoll clas- 
sification scheme have recently been e xtensively studied 
and discussed by ISchnerr et alJ (|2003h . They found that 
its motion through the HID is in the opposite sense to 
that in the CD, with the X-ray count rate increasing again 
when the source moves into the spectrally hard state, con- 
trar y to most atoll and Z sources (see our Fig. F2] and 
also Wiinands et al. (1997) for similar behavior in the Z 



source Cyg X-2). Moreover, the strength of the very low 
frequency variability also changes in the opposite sense to 
that observed in Z and atoll sources. 

In an attempt to fit GX 13+1 within th e Z/atoll 
scheme as an atoll source, Sc hnerr et al ] ll2003l) tried to 
explain the source's peculiar behavior with the presence 
of a relativistic jet with an axis nearly aligned to our line of 
sight. The radio emission could then be boosted by more 
than an order of magnitude, putting the radio emission 
of GX 13+1 in accordance with measurement and upper 
limits of other atoll source; the unusual X-ray phenomena 
could be the consequence of a better view of the X-ray 
emitting regions associated with the base of the jet, pos- 
sibly assisted by Doppler boosting. 

As we showed above, the X-ray and radio luminosi- 
ties of GX 13+1 are actually in the range expected for 
Z sources, without requiring any unusual jet geometries, 



making a Z source nature more likely in our opinion (the 
atoll sources GX 9+1, GX 9+9 and GX 3+1 have also 
similar X-ray luminosities but much lower radio luminosi- 
ties). This is supported by our spectral fit parameters, 
which a re quite similar to those found in the Z sources 
(see e.g. ldi Salvo et all2000HDi Salvo et al.l200ll) and less 
like those in the atoll sources (see e.g.lDi Salvo et al . 2000; 
iBarret et alJl2000t iBarret fc Olivell2002|) . and by the fact 
that GX 13+1 rarely shows type I X-ray bursts. Based 
on the morphology in our CD GX 13+1 was probably ob- 
served in the normal branch and flaring branch. Although 
no clear indications for normal branch o scillations are 
found , the band limited noise measured bv lSchnerr et all 
( 2003) has properties s imilar to the norma l branch oscilla- 
tions found in GX 5 -1 ijjonker et alJl2002h and GX 340+0 
l(,Tonker et alJ l200nl> close to the normal branch/flaring 
branch vertex, where it is very broad (Q < 1) and rather 
weak (<2% rms). In GX 5-1 and GX 340+0 this broad 
bump evolves into a strong 6 Hz QPO as the spectrum 
hardens - this is not observed in GX 13+1. There are ad- 
ditional differences with the other Z sources, in particular 
with respect to the behavior of the low frequency variabil- 
ity. 

Finally, we note that although the source was more 
similar to the Z sources than to the atoll sources during 
our observations (at a luminosity three times as high as 
during the EXOSAT era), this does not exclude that at 
lower luminosities it behaves more atoll-like. 



5. Summary 

Our simultaneous X-ray/radio observations of GX 13+1 
revealed a strong dependence of the radio brightness on 
the X-ray state of the source. In the hard spectral state, 
which was at least 4 times brighter in the radio than the 
soft state, we also found a correlation between the X-ray 
and radio properties on a short time scale, with changes 
in radio having a delay of ~40 minutes with respect to 
those in X-rays. We attribute this delay to the time it 
takes for the flow to become optically thin in the radio. 
The absence of strong dips in the X-ray light curve during 
the radio flare suggests that only a small amount of the 
matter is redirected from the inflow to the outflow. On the 
basis of a comparison with atoll and Z sources we conclude 
that the source is more similar to the Z sources, although 
several properties of GX 13+1 remain unexplained. 
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